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Abstract 

' We construct the covariant nonlocal action for recently suggested 



long-distance modifications of gravity theory motivated by the cos- 
mological constant and cosmological acceleration problems. This con- 



0*N struction is based on the special nonlocal form of the Einstein-Hilbert 

action explicitly revealing the fact that this action within the covariant 
<^> curvature expansion begins with curvature-squared terms. 

New approach to the solution of the cosmological constant problem con- 
On! sists in the assumption that, instead of adjusting the vacuum energy of 
r-| ! quantum matter to zero (or to small value of the cosmological acceleration), 
one should modify the purely gravitational sector of the theory in far in- 
frared region. Matter sources with wavelengths comparable with the horizon 
size of the present Universe L ~ 1/Hq ~ 10 28 cm gravitate with the long- 
distance gravitational constant Gld which is much smaller than the conven- 
tional Planckian value Gp. Therefore, the vacuum energy £, T^ u = Sg^, 
of TeV or even Planckian scale (necessarily arising in all conceivable models 
with spontaneously broken SUSY or in quantum gravity) will not generate a 
catastrophically big spacetime curvature incompatible with the tiny observ- 
able Hubble constant Hq ~ Gld& *C GpS. This mechanism is drastically 
different from the old suggestions of supersymmetric cancellation of £ [1], 
because it relies on the fact that the nearly homogeneous vacuum energy 
gravitates very little, rather than it is itself very small. 

This idea implies that the gravitational coupling constant should be pro- 
moted to the level of the operator, Gp => G(0), which for sake of covariance 
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can be regarded as a function of the covariant d'Alembertian □ = g al3 V a 'V p, 
interpolating between the Planck scale of the gravitational coupling constant 
G P = 1/lQirMp for local matter sources of size <C L and the long distance 
gravitational constant Gld — G(0) with which the sources nearly homoge- 
neous at the horizon scale are gravitating [2]. The modified equations of 
motion were suggested to have the form of Einstein equations 

m\u)(r^- 1 - 9 ^ =hr m (l) 

with "nonlocal" inverse gravitational constant or Planck mass, 

' M 2 (D) = Mp (1 + JF(L 2 D )), (2) 



16ttG(D) 



being some function of the dimensionless combination of □ and the additional 
scale L - the length at which infrared modification becomes important. If 
the function of z = L 2 satisfies the conditions, T{z) — > at z 3> 1, and 
T{z) — > .F(O) ^> 1 at z — > 0, then the long-distance modification is inessential 
for processes varying in spacetime faster than 1/L and is large for slower 
phenomena at wavelengthes ~ L and longer. 

One difficulty with this construction is that for any nontrivial form factor 
JF(L 2 D ) the left hand side of (1) does not satisfy the Bianchi identity and, 
therefore, cannot be generated by generally covariant action. Obviously, this 
makes the situation unsatisfactory because of a missing off-shell extension 
of the theory, problems with its quantization, etc. Here we suggest to cir- 
cumvent this problem by resorting to the weak field approximation, which is 
certainly valid in the infrared regime. This means that Eq. (1) should be 
understood only as a first term of the perturbation expansion in powers of 
the curvature. Its left hand side should be modified by higher than linear 
terms in the curvature, and the modified nonlocal action Snl[9] should be 
found from the variational equation 

= M P g^(l + T^U )) (rT - \<TR) + O [RU- (3) 

Flexibility in higher orders of the curvature allows one to guarantee the 
integrability of this equation and to construct the nonlocal action as a gener- 
ally covariant (but nonlocal) curvature expansion. Here we explicitly present 
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this construction along the lines of covariant curvature expansion of [3]. As 
a starting point we consider a special nonlocal form of the Einstein-Hilbert 
action revealing its basic property - the absence of a linear in metric pertur- 
bation part (on flat-space background). Then we introduce a needed long- 
distance modification by a simple replacement of the nonlocal form factor in 
the curvature-squared term of the obtained action [4]. The paper is accom- 
plished by a discussion of the nature of nonlocalities in quantum-gravitational 
and brane- induced models of [5]. In particular, the fact that curvature ex- 
pansion for the action begins with the quadratic order is revisited from the 
viewpoint of the running gravitational coupling constant and the issue of 
acausality of nonlocal effective equations, raised in [2], is reconsidered. 
The Einstein action in the Euclidean asymptotically-flat spacetime 



S E [g] = -M 2 P / dx g 1 ' 2 R(g)+M 2 P da» {&> V - d^h) . (4) 




includes the bulk integral of the scalar curvature and the Gibbons-Hawking 
surface integral over spacetime infinity, \x\ — > oo. The latter is usually called 
the Gibbons-Hawking action which in the covariant form contains the trace 
of the extrinsic curvature of the boundary (with an appropriate subtraction 
of the flat space background contribution). This surface term guarantees the 
consistency of the variational problem for this action which yields as a metric 
variational derivative the Einstein tensor. 

The action (4) is explicitly linear in the curvature, but this linearity is 
misleading, because its variational derivative — the Einstein tensor — is also 
linear in the curvature. Therefore, it is at least linear in metric perturbation 
on flat-space background, ~ h^ u , and the perturbation theory for Se [g] 
should start with the quadratic order, O [ ] ~ O [ R 2 ] . This is a well 
known fact from the theory of free massless spin-2 field. Our goal is to make 
this /^-expansion manifestly covariant by converting it to the covariant 
curvature expansion. A systematic way to do this is to use the technique of 
covariant perturbation theory of [3], which begins with the derivation of the 
expression for the metric perturbation in terms of the curvature. 

Expand the Ricci curvature in h„ v on flat-space background 






1 □ v + 1 {y,x v + v„x„) + o [ h\ v ], 



(5) 



Vfj,h, and solve it by iterations as a nonlocal expansion in 
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powers of the curvature. This expansion starts with the following terms 

v = + + v„/„ + o (6) 

Here 1/D denotes the action of the Green's function of the covariant metric- 
dependent d'Alembertian on the space of symmetric second-rank tensors with 
zero boundary conditions at infinity and the term V fl f u + V u f^ in (6) reflects 
the gauge ambiguity in this solution (it originates from the harmonic-gauge 
terms in the right-hand side of (5)). 

Now restrict ourselves with the approximation quadratic in R^ v (or equiv- 
alently h^) and integrate the variational equation for Se [<?]■ Since the vari- 
ational derivative is at least linear in h^ v , 5SE/8g^ v ~ h a/3 , the quadratic 
part of the action in view of this equation is given by the integral 

S E [g] = \f dx VW^j + O [RU- (7) 

Substituting the Einstein tensor for SSe/Sq^ and (6) for and integrating 
by parts one finds that the contribution of the gauge parameters f ^ vanishes 
in view of the Bianchi identity for the Einstein tensor, and the final result 
reads 

S E [g] = MlJ dxg 1 ' 2 { -{RT -\g^R) ^ + O [R^] } . (8) 

This is the covariant nonlocal form of the local Einstein action [6, 7, 4]. 
This nonlocal incarnation of (4) explicitly features the absence of the linear 
in curvature term, which can be clarified by the subtraction effect of the 
Gibbons-Hawking term. 

In asymptotically-flat (Euclidean) spacetime with the asymptotic behav- 
ior of the metric g^ = S^ + h^, = O (l/|:c| d_2 ) , \x\ — > oo, the Gibbons- 
Hawking term in Cartesian coordinates can be transformed to the bulk inte- 
gral of the integrand d lx {d v h ilv — d^hj - the linear in h^ v part of the scalar 
curvature. Similarly to the above procedure this integral can be covariantly 
expanded in powers of the curvature. Up to quadratic terms inclusive this 
expansion reads [4] 

Jja^h^-d,h)= j dxg 1 / 2 ^R-(R^ - l -g^R)^R^ + ..^. (9) 
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As we see, when substituting to (4) its linear term cancels the Ricci scalar 
part and the quadratic terms reproduce those of (8). Obviously, this type of 
expansion can be extended to arbitrary order in curvature. 

Long distance modification of the Einstein action that would generate 
(3) as the left-hand side of the gravitational equations of motion now can 
be simply obtained from the nonlocal form of the Einstein action (8). It 
is just enough to make the following replacement in the quadratic part of 
(8), 1/D — > (1 + JF(L 2 D))/D. Indeed, the variation of the Ricci tensor here 
and integration by parts " cancel" □ in the denominator. All commutators of 
covariant derivatives with the □ in JF(L 2 D) give rise to the curvature-squared 
order which is beyond our control. This recovers the Einstein tensor term of 
(3) with the needed "nonlocal" Planckian mass Mp(l + JF(L 2 D )). 

Thus, the long-distance modification in question takes the form 

Snl[ 9 ] = -J dxg 1 ' 2 {(RT - \g^R)^^R, v + 0[R%}} . (10) 

It is manifestly generally covariant, and its variational derivative (the left 
hand side of the modified Einstein equations) exactly satisfies the Bianchi 
identity and does not suffer from the concerns of [2]. This action is not 
unique though, because it is defined by a given form factor JF(L 2 D ) only in 
quadratic order, while we do not have good principles to fix its higher-order 
terms thus far. 

One of the main mechanisms for nonlocalities of the above type is the 
contribution of graviton and matter loops to the quantum effective action. 
In quantum theory the concept of a nonlocal form factor replacing a cou- 
pling constant is not new. In fact this concept underlies the notion of the 
running coupling constants and sheds new light on the cosmological constant 
problem also from the viewpoint of the renormalization theory. For sim- 
plicity, consider QED or Yang-Mills theory in the quadratic order in gauge 
field strength F 2 V . The transition from classical to quantum effective action, 
S — > S e s, boils down to the replacement of the local invariant by 

g- 2 J dxFl -> J dxF„ v g;£{-n)F^. (11) 

Here the effective coupling constant g^{— □) is a nonlocal form factor which 
can be obtained from the corresponding solution of renormalization-group 
equation. Obviously, it plays the role of T{L 2 U) above. 
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This concept, however, fails when applied to the gravitational theory in 
the sector of the cosmological and Einstein-Hilbert terms 1 . Indeed, naive 
replacement of ultralocal cosmological and gravitational coupling constants 
by nonlocal form factors, 

J dx g 1/2 (A - Mp R) -> J dx g 1/2 (k{U) - M 2 P (n) R) , (12) 

is meaningless because the action of the covariant d'Alembertian on the right 
hand side always picks up its zero mode, and both form factors reduce to 
their numerical values in far infrared, A(0), Mj>(0). Therefore, even if one 
has solutions of renormalization group equations for A and Mp, like those 
of [8] , one cannot automatically recover the corresponding pieces of effective 
action or the corresponding nonlocal correlation functions. 

The construction above suggests that the running coupling constant "de- 
localization" of Mp should be done in the (already nonlocal) representa- 
tion of the Einstein action (8). As its curvature expansion begins with the 
quadratic order, one can insert the nonlocal form factor Mj>(D) between two 
curvatures so that no integration by parts would result in its degeneration 
to a trivial constant. It would be interesting to see how a similar mechanism 
works for the nonlocal cosmological "constant" A(D). 2 The mechanisms of 
its generation due to infrared asymptotics of the effective action, or late-time 
asymptotics of the corresponding heat kernel, are discussed in [12]. 

Nonlocalities of the type (10) also arise in a certain class of braneworld 
models [13, 5]. They cannot appear in models of the Randall-Sundrum 
type with strictly localized zero modes, because in these models nontriv- 
ial form factors basically arise in the transverse-traceless sector of the action 
(as kernels of nonlocal quadratic forms in Weyl tensor [7]). In contrast to 
these models, the nonlocal part of (10) is not quadratic in the Weyl ten- 
sor, / dx g l l 2 W 2 yafj ~ J dxg l / 2 (R 2 lv — \R 2 ) (with the insertion of a nonlocal 
form factor between the curvatures). Rather, (10) includes the structure 

1 When applied to formally renormalizable (albeit non-unitary) curvature-squared grav- 
itational models [8]. 

2 One should expect that a quadratic action for the cosmological term would read as 
A J dxg 1 l 2 R! lv {\/n 2 )R fll ,. This structure (also suggested in [9] and discussed within the 
renormalization group theory) appears in two-brane models [7] and expected as a covariant 
completion of the mass term in models of massive gravitons [10] and discussions of the 
van Damm-Veltman-Zakharov discontinuity [11]. 
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/ dx g 1 ^ 2 {R 2 lu — \R 2 ) which contains the conformal sector. It is this sector 
which is responsible for the potential resolution of the cosmological constant 
problem. It becomes dynamical in models with metastable graviton like the 
Gregory-Rubakov-Sibiryakov model [13] or brane induced gravity models of 
the Dvali-Gabadaze-Porrati (DGP) type [5]. In particular, for the (4+1)- 
dimensional DGP model with the bulk Gab{X) and induced on the brane 
g^{x) metrics 

S DGP [G] =M 3 J d^XG 1 ' 2 b R{G AB ) + M 2 P J d'xg 1 ' 2 A R{g, v ) (13) 

(we disregard here relevant Gibbons-Hawking terms) the effective braneworld 
action takes the form (10) with the form factor JF(L 2 D) which is singular at 
□ -> [5, 14] 

^ = Z^B' £ = §- < 14 > 

where M ~ 10~ 21 Mp ~ 100 MeV is a mass scale of the bulk gravity as 
opposed to the Planckian scale of the Einstein term on the brane Mp ~ 10 19 
GeV. This model is interesting, because its approximate cosmological FRW 
equation of motion, 



H 2 -- = ^, (15) 
L 6Mp' V ; 



admits the self-accelerating regime with H — 1/L at late stages of evolution 
with matter density p — > [14]. Unfortunately, this model can hardly serve as 
a consistent candidate for the description of our Universe because of the pres- 
ence of the sufficiently low strong-coupling scale (Mp/L 2 ) 1 / 3 ~ 10~ 8 cm _1 , 
which invalidates its applications to phenomena already at distances smaller 
than 1000 km [15] (see, however, Ref. [16] advocating that this difficulty can 
be circumvented by special type of regularization). 

Form factors like (14) are unambiguously defined only in the Euclidean 
space with the negative semi-definite d'Alembertian □. This raises the prob- 
lem of their continuation to the Lorentzian spacetime where the issues of 
causality and unitarity become important. The principles of this continuation 
depend on the physical origin of nonlocality in JF(L 2 D). In particular, scat- 
tering problems for in-out matrix elements, ( out \ip\in), in spacetime with 
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asymptotically-flat past and future imply a usual Wick rotation. The prob- 
lem for in- in mean value of the quantum field, (f) — ( in \ tp \ in ), is more com- 
plicated and incorporates the Schwinger-Keldysh diagrammatic technique 
[17]. In this technique the effective equations for <fi cannot be obtained as 
variational derivatives of some one-field action functional. However, there 
exists a special case of the initial quantum state — the Poincare-invariant 
in-vacuum in asymptotic past, \in) = \ in, vac). Effective equations for 
in this vacuum can be obtained by the following procedure [3]. Calcu- 
late the Euclidean effective action in asymptotically-flat spacetime, take its 
variational derivative containing the nonlocal form factors which are uniquely 
specified by zero boundary conditions at Euclidean infinity. Then formally go 
over to the Lorentzian spacetime signature with the retardation prescription 
for all nonlocal form factors. These retarded boundary conditions uniquely 
specify the nonlocal effective equations and guarantee their causality. This 
procedure was proven in [3] and also put forward in a recent paper [18] as 
the basis of the covariant nonlocal model of MOND theory. 

This procedure justifies the Euclidean setup used above and suggests 
that in this setting no contradiction arises between the nonlocal nature of 
the Euclidean action and causal nature of nonlocal equations of motion in 
Lorentzian spacetime. In this respect the situation is essentially different 
from the assumptions of [2] where acausality of equations of motion is nec- 
essarily attributed to the nonlocal action. No such assumptions are needed 
in effective equations for expectation values which are fundamentally causal 
despite their nonlocality. These equations have interesting applications in 
quantum gravitational context and, in particular, show the phenomenon of 
the cosmological acceleration due to infrared back-reaction mechanisms [19]. 

The situation with brane induced nonlocalities and their causality status 
is more questionable and conceptually open. For example, the branch point 
of the square root in the nonlocal form factor (14) is apparently related to 
different branches of cosmological solutions of (15) including the scenario of 
cosmological acceleration [14]. Therefore, in contrast to tentative models 
of [2] with finite .F(O) ^> 1, which only interpolate between two Einstein 
theories with different gravitational constants Gld ~ Gp/JF(0) <C Gp, the 
DGP model is anticipated to suggest the mechanism of the cosmological ac- 
celeration. This implies the replacement of the asymptotically-flat spacetime 
by the asymptotically-deSitter one. For small values of asymptotic curva- 
ture (as is the case of the observable horizon scale Hl/Mp ~ 10~ 120 ) the 
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curvature expansion used above seems plausible, although the effect of the 
asymptotic curvature might be in essence nonperturbative. Therefore, the 
above construction might have to be modified accordingly. In particular, the 
expansion in powers of the curvature should be replaced by the expansion 
in powers of its deviation from the asymptotic value — hg^Roo (d is 
a spacetime dimensionality). This would introduce in the formalism as a 
free parameter the value of the curvature in far future, Roo ~ Hq, reflecting 
the measure of acausality in the model. The resulting modifications in the 
above construction are currently under study and will be presented elsewhere. 
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